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We investigate the properties of crystalline phantom membranes, at the crumpling transition and 
in the flat phase, using a nonperturbative renormalization group approach. We avoid a derivative 
expansion of the effective average action and instead analyse the full momentum dependence of the 
elastic coupling functions. This leads to a more accurate determination of the critical exponents and 
further yields the full momentum dependence of the correlation functions of the in-plane and out- 
of-plane fluctuation. The flow equations are solved numerically for D = 2 dimensional membranes 
embedded in a d = 3 dimensional space. Within our approach we find a crumpling transition 
of second order which is characterized by an anomalous exponent rjc ~ 0.63(8) and the thermal 
exponent f ~ 0.69. Near the crumpling transition the order parameter of the flat phase vanishes 
with a critical exponent j3 « 0.22. The flat phase anomalous dimension is 77/ « 0.85 and the 
Poisson's ratio inside the flat phase is found to be af ~ —1/3. At the crumpling transition we flnd 
a much larger negative value of the Poisson's ratio ctc « —0.71(5). We discuss further in detail the 
different regimes of the momentum dependent fluctuations, both in the flat phase and in the vicinity 
of the crumpling transition, and extract the crossover momentum scales which separate them. 

PACS numbers: 05.20.-y,11.10.Lm,68.65.Pq,05.10.Cc 



I. INTRODUCTION 



Crystalline membranes have a fixed connectivity of 
their constituent particles and therefore a finite in-plane 
shear modulus which distinguishes them from fluid mem- 
branes [H]. Physical crystalline membranes are two di- 
mensional membranes embedded in three dimensional 
space and the embedding allows for out-of-plane fluctu- 
ations into the third dimension which are absent in a 
truly two-dimensional crystal. The coupling of the in- 
plane modes to the out-of-plane ones is responsible for 
a strong renormalization of the infrared (IR) character 
of both modes. Both acquire large anomalous dimen- 
sions which render the fluctuation of the local normals 
of the membrane finite and thus stabilize a flat phase 
in which the normals are ordered 0. The physics of 
such flat crystalline membranes has recently received re- 
newed attention because of the discovery of graphene 
3, the thinnest possible crystalline membrane Q. At 
higher temperatures, the flat phase eventually becomes 
unstable and looses its orientatiqnal order of the nor- 
mals at the crumpling transition |5|. Several approaches 
have been used to investigate the crumpling transition 
of phantom membranes. In contrast to real membranes, 
phantom membranes have no self-avoiding contact inter- 
action which makes them more accessible to analytical 
approaches. Paczuski et al. studied a generalized model 
of a _D-dimensional elastic manifold embedded in a d- 
dimensional space within a leading order e = 4 — D ex- 
pansion and found for I? = 2 (e = 2) a first order tran- 
sition for all d < 219. The accuracy of this approach, 
which is well controlled only for small e, at e = 2 is 
however in doubt. The results from the e expansion dis- 
agree for Z? = 2 with the results of a large d expansion 



0| and its extension, the selfconsistent screening ap- 
proximation (SCSA) |8l-[l0| where a continuous transition 
was found for D = 2 and d = 3. Further support for a 
continuous transition came from an elastic model with a 
constraint arising from infinitely large coupling constants 
for the stretching modes [ll|, from a Monte Carlo (MC) 
renormalization group analy sis fl2| . and from early MC 
simulations, see references in [l^. More recent MC simu- 
lations found however evidence for the coexistence of two 
separate phases at the critical temperature which would 
be consistent with a first order transition [3 113; see 
also [l6l . [l7| . A recent nonperturbative renormalization 
group (NPRG) analysis, based on a derivative expan- 
sion of the effective average action, could reproduce both 
the leading order results of the e = 4 — D expansion as 
well as the leading order of a large d expansion, within a 
unified framework [l8|. The technique has recently also 
been applied to anisotropic membranes [l^. For D = 2 
and d = 3 evidence of a second order phase transition 
was found, but a weak first order transition could not be 
ruled out because of a weak dependence of the flow on 
the employed cutoff scheme. 

Crystalline membranes are characterized by an anoma- 
lous elasticity in which none of the local elastic constants 
remain finite. The elasticity becomes fully nonlocal in 
the IR limit and the elastic coupling constants become 
coupling functions with a momentum dependence which 
for asymptotically small momenta is characterized by 
anomalous exponents. Here we present a NPRG anal- 
ysis where the full momentum dependence of the cou- 
pling functions is kept, which allows for a more accurate 
analysis of the crumpling transition which we find to be 
continous. Our approach further allows to calculate the 
thermal fluctuations of the membrane beyond the asymp- 
totic regime. We recover the results of Ref. [31 if the 
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momentum dependence of the coupling functions is ne- 
glected. We previously applied our approach to compute 
the thermal fluctuations of free standing graphene [l^l 
and found excellent agreement for all momenta with MC 
simulations i21j, 22^ . This analysis allowed to identify the 
characteristic scale of ripples in free standing graphene Q 
with the Ginzburg scale of the nonlinear elastic theory of 
crystalline membranes. The SCSA is also capable of ac- 
cessing fluctuations at finite momenta and was applied to 
investigate the thermal fluctuations of graphene [10], but 
proved to be, at finite momenta, less accurate than the 
NPRG approach. Here, we extend our approach and in- 
vestigate in detail the behavior of a crystalline membrane 
inside the flat phase and at the crumpling transition crit- 
ical point. 



All approximations of our approach are included in 
our nonlocal ansatz of the effective average action and 
not at the level of the flow equations for irreducible ver- 
tices. Once the effective average action, which respects 
the full symmetry of the symmetric phase, is specified, 
the flow equations of the irreducible vertices are uniquely 
determined and can be solved exactly. In all previous 
schemes to calculate the full momentum de pen dence of 
self-energies of classical models (see Refs. |2j] for ap- 
proaches for the $^ theory), some approximations were 
done in addition to the truncation of the effective average 
action which in principle can lead to a violation of Ward 
identities. Our approach is quite general and could also 
be applied to other physical models. Similar flow equa- 
tions were derived for a model of interacting bosons [24| . 
However, additional approximations had to be employed 
to solve the flow equations in which case the approach be- 
comes equivalent to that presented in [l^ . The approach 
presented here assumes an ordered state but allows to ap- 
proach and analyse the critical point from within the low 
temperature phase. 



The structure of this manuscript is the following. 
In Sec. In] we introduce and define the usual Landau- 
Ginzburg model of crystalline phantom membranes, 
which are membranes with no self-avoiding contact inter- 
action. In Sec, mil we present the derivation of our NPRG 
flow equations for the nonlocal elastic coupling functions. 
Expressions for the two longest diagrams of the flow equa- 
tions can be found in the appendix. In Sec. lIVI we present 
results from a completely self-consistent numerical solu- 
tion of the flow equations. We analyse the different scal- 
ing regimes in the momentum dependence of the out-of- 
plane and in-plane fluctuations and determine the critical 
exponents of the membrane both in the flat phase and at 
the crumpling transition fixed point. We further discuss 
the Poisson's ratio, which is different at the crumpling 
transition and inside the flat phase. Finally, in Sec.lVlwe 
present a summary and conclusion of this work. 



II. LANDAU-GINZBURG MODEL 

Our starting point is the usual Landau- Ginzburg 
model for crystalline phantom membranes (which have 
no contact interaction) @, with H = + The 
bending part of the membranes is described by (silent 
indices are summed over) 



d^'xidadaRY 



(la) 



and the stretching is described by 



{daRf 



{daR ■ dbR) 



--{daR-daRf] , 



(lb) 



where i? is a D -f 1 dimensional vector parameterizing 
the D dimensional membrane which is embedded in a 
d — D + \ dimensional space. For notational conve- 
nience we restrict our analysis to c? = D -I- 1 , it is however 
straightforward to extend our analysis to arbitrary d. In 
terms of the derivatives rUa — daR the stretching part 
becomes 



n"^ = / d^x 



ro 



• rria + ^("^a ' rUb)^ 



+ - (ma ■ ma) 



(2) 



and one can already anticipate a transition near f o c± 
from a flat configuration, which exists for negative and 
where the fields ma acquire a finite expectation value 
"ma.o — (rria) 7^ 0, to a crumpled phase with (ma) — 0. 
The stretching part Eq. ([2]) is appropriate for an isotropic 
membrane whose consituents particles have a fixed con- 
nectivity (such membranes are also refered to as tethered 
or polymerized membranes). A crystalline membrane has 
in general a lower symmetry and additional terms quartic 
in ma are allowed, except for a two-dimensional hexago- 
nal lattice 6] . We here concentrate on models described 
by Eq. appropriate for isotropic manifolds or mem- 
branes with hexagonal lattices and D — 2. 

For the flat phase, we introduce J 7^ as the magni- 
tude of the order parameter which is defined as {R{x)') = 
JxaBa with o = 1 . . . D and are orthonormal vec- 
tors which span the plane of the membrane. Defining 
Uab = {rUa ■ mb - maS) ■ rnbs))l2 and maX) = Jea, one 
can write [6| 



d'^x 



2 aa 
t2i~ 



+ j d^a;(fo + J'[n + D\l2]'jUaa ■ (3) 

The last term has to vanish to guarantee thermodynamic 
stability and we thus have the mean field result 



•J \ 



(4) 
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where the subscript Aq indicates that this is a relation 
among the bare parameters which are assumed to be de- 
fined at some scale Aq. 



III. DERIVATION OF THE 
NONPERTURBATIVE RENORMALIZATION 
GROUP FLOW EQUATIONS 

The NPRG approach is based on the exact flow equa- 
tion of the effective average action Fa [1^ 



dTA 1 



Ra 



dA 



(5) 



where the fields </>, (/>' here represent the components of 
the vector field R and the trace stands for a momentum 
integral and a sum over internal indices. The function 
is a regulator that removes IR divergences arising from 
modes with A; < A and will be specified below. 

In principle, since the NPRG is not based on an expan- 
sion in a small parameter, the NPRG is not a-priori con- 
trolled. It has however been applied with much success 
to critical phenomena [27| and it usually reproduces lead- 
ing order results of pcrturbative RG techniques such as 
e expansions around the upper (or lower) critical dimen- 
sion. For the case of crystalline membranes, the simplest 
possible NPRG ansatz [l^ already reproduces the lead- 
ing order results of both an e = 4 — D expansion as well 
as a 1/d expansion. The approach presented here allows 
to calculate the momentum dependence of the membrane 
fluctuations which were found to be in perfect agreement 
with large scale MC simulations of unsuspended graphene 
(20j . The NPRG technique thus seems to be quite reliable 
when applied to crystalline membranes. 



A. Nonlocal elasticity 

Our approximation of the effective average action con- 
sists of a nonlocal bending part F*^, which is quadratic in 
the fields, and a nonlinear stretching part F^*, which also 
includes quartic terms. 



L 
2 



i J d^x d^x'kAix - x')dlR{x)dlR{x') 

flA{x - x')Uab{x)Uab{x') 



d^x d^ X 



Ir 



+ -\a{x - x')Uaa{x)Ubb{x') 



(6) 



This ansatz is simply a nonlocal generalization of the bare 
model defined via Eqs. ([Ta]) and Nonlocal correla- 
tions are known to be dominant in the IR limit of crys- 
talline membranes where the Fourier-transformed cou- 
pling functions scale anomalously. Both in the flat phase 
and at the critical point of the crumpling transition, one 



has a pronounced hardening of the out-of plane fluctua- 
tions and a softening of the in-plane fluctuations. This 
is expressed by a divergent form of the bending coupling 
function and a vanishing of the bulk and shear modulus. 
In the flat phase one has k{q) ~ q~^f whereas the elastic 
coupling functions vanish as jl{q) ^ Ml) ^ <f ■ Here, 
•qf and 77" are the anomalous exponents of the out-of 
plane fluctuations and the in-plane fluctuations, respec- 
tively. Note that these anomalous exponents are in fact 
not independent. Invariance of the original model under 
rotations implies a relation (Ward identity) among the 
anomalous dimensions [dl, 



4 - D - 2r?/ . 



(7) 



Since all fluctuations become anomalous in the IR limit, 
neither K(q), /i(fc), nor A(fc) are analytic functions in the 
limit A — i> and all fluctuations become nonlocal, a sit- 
uation familiar from the behavior at a critical point of a 
continuous phase transition. In a leading order derivative 
expansion only the momentum independent parts of the 
coupling functions are kept, which is sufficient to extract 
their asymptotic momentum dependence by approximat- 
ing fA=o{k) ~ /A=fe(0) where / is any of the three cou- 
pling functions. We choose a more accurate approach, 
which also captures the corrections to the asymptotic be- 
havior, and simply keep the full momentum dependence 
of the coupling functions for all A. 

It is convenient to write the flow equation ([5]) of the 
functional Fa as a flow equation of vertex functions, 
which can be obtained from a field expansion of both 
sides of ([5]). Since we are interested in the symme- 
try broken phase, we work with fields corresponding to 
the in-plane fluctuations u and out-of-plane fluctuations 
h rather than the original R fields. These are intro- 
duced via the fluctuations of rUa around the ordered 
state, Ama = daR ~ ruafi with Ama = {daU,dah) and 
Tna = Ja^u- We therefore write the expansion 



TA[u,h] = 



E 



1 



7, ' rr?, ' 



T 



(n+m) 
A.,h . . . hai. 



(qi,...,q„;pi,...,p„) 



X h, 
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h 7/"i 



where the momentum integrals are defined as 



(8) 



(9) 



The irreducible vertices entering Eq. ([SJ can be re- 
lated to correlation functions. The simplest example is 
the Dyson equation which relates the one particle Green's 
function to the irreducible self-energies. The Dyson equa- 
tion for the Green's function Ghh of the h field, with 



(hqh^q,) = VSq,q'Ghh{q) , 



(10) 
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where V is the D-dimensional volume, is 

Gr,l{q)=G,XiQ) + ^hHiq)- (11) 

The self-energy is of the form (here and below we sup- 
press in our notation the explicit A dependence of the 
coupling parameters kg flq, and Xq) 



^hhiq) = rf^^(q, -q) = (Kg - KAjg'' 



(12) 



and the cutoff dependent non-interacting Green's func- 
tion is 



GoXil) = HaoQ^ + RA{q), 



(13) 



where denotes the bare and momentum independent 
value of the initial coupling constant k defined at the 
ultraviolet (UV) cutoff Aq. Here, we introduced the reg- 
ulator function R/^{q), which regulates the IR limit of 
the propagator in such a way that the IR divergence at 
g — ^ is removed at finite A. Since we will later solve the 
NPRG flow equations numerically, we choose an analytic 

cutoff fii, 



q 



exp[(g/A)4] - 1 



where 



~(o) _ ~ 



(14) 



(15) 



is the q = component of the Fourier transform of k\(x). 

For finite A, we then have limq_>o = itj^' A"* and Go. a 
is non-divergent for q ^ 0. 

The Green's functions of the in-plane modes, defined 
via 



(16) 



can be written in terms of transverse (_L) and longitudinal 
(II) components using the projectors 

P,i{k)^Sab-k''k''/k\ (17a) 

Pi(fe) = k^k'/k' . (17b) 

This yields 

Gatik) = G^(fc)Pi(fe) + G[|(fc)Pi(fc) (18) 

with G~^ = Gq \ + Sq for a =±, || and where the self- 
energies are defined as projections of 



Safc(fc)-rf,^(fc,-fe) 

With the two-point irreducible vertex 



(19) 



r^',)(fci, -fca) = {Sabki ■ k2 + k^k\) -f Xk^k^kl^ 



(20) 



one finds the in-plane self-energies 

E_l(A:) = Jlflkk'^ + {kk - fiA„)k'^ , 



(21a) 



= Jli2fik + ~Xk)k^ + [kk - kA,)k\ (21b) 

The effective average action Eq. ^ contains further 
three- and four-point vertices of the form (we use k for 
momenta of w-fields and p for momenta of /i-fields) 



rifa.a3(fcl.fe2,fe3) 



-iJa 



J2 {2flkAkc. ■ kp)k^'' 



+ XkAkp-k^)k^''}da,^a, (22a) 



ril(pi,P2; k) = -iJA{Afc [iPi ■ k)p^2 + iP2 ■ fc)K] 



Afc(pi-P2)r} 



(22b) 



rit^...a.(fcl---fc4) = ^ E {^fik^A^a ■ k,){kp ■ ks) 

— "^1234 

+ Afc„^ (fe„ • kp){k^ ■ ks)}Sa^,ai3Sa^,as : (22c) 



rLaia2(Pl'P2;fel,fe2) = { [(Pl " ^1 ) (P2 ' ^2) 
+ (Pl • fe2)(P2 • ki)] + Ap,,(pi ■P2)(fcl • k2)}6a,.a2 : 

(22d) 



= 8 E {2/ip„^ + V^}(Pa •P7)(P/3 -Ps): 

(a/97(5)=7'i234 

(22e) 

where k^p = I^q + k^l and J2(afi'y)=Vi23 denotes a sum- 
mation over all permutations of (123). The subscript h 
refers to /i-fields while subscripts a and refer to u- 
fields. Since we neglect all irreducible correlations not 
explicitely included in Eq. higher order vertices, with 
more than four legs, do not appear in our approximation. 
Finally, we shall also need the single scale propagators 
which are defined via 



Gab = -(GiPi, + G|Pi)aAi?A. 



Ghh — ^G^,j9a-Ra ■ 



(23a) 
(23b) 



B. Flow equations for momentum dependent 
vertices 

We will now derive the flow equations for the order 
parameter Ja and the coupling functions Kg, flq and Xq. 
The flow of can be extracted from the flow equation 
of the one-point vertex of the M-flelds whereas the other 
three coupling functions can be extracted from the flow 
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equations of the self- energies. All higher order vertices 
presented in Sec. lIII Al are completely determined by J|, 
Kg, jlq and Xq, and the flow equations for the one- and 
two-point vertices are thus closed. 

The only approximation of our approach is the ap- 
proximation for the effective average action as expressed 
through Eq. (|6]) and in the derivation below no further 
approximation is required. The resulting flow equations 
are uniquely determined by Eqs. ([5]) and (|6]) and obey 
the full symmetry of Eq. © . 

The NPRG flow of the one-point function ri^-* is 



1 



.(3) 



(24) 



where we used the fc-space representation of the order 
parameter field 



d 



(25) 



'\r\j^\j\Z = - 1/2 - 1/2 



FIG. 1: (Color online) Diagrammatic representation of the 
flow equation which determines the flow of the order parame- 
ter and which arises from the vanishing flow of the one-point 
vertex [2^. Wavy lines correspond to Gab propagators and 
solid lines to Ghh- Lines with a dash correspond to single- 
scale propagators. The small open circle with a dot denotes 
a derivative of the order parameter with respect to A. 

The condition dA^'a\k) — then yields an equa- 
tion which is shown diagrammatically in Fig. [1] and from 
which the flow of the order parameter can be determined. 
One finds 



Kt 



J dq q^+^[{D-l)G^{q) 



[2/io + DXo]D 

X [2fig + 2/io + DXa] 

+ G|i {q) \2{2fig + fio) + 2Xq + DXo 



Ghhiq)[2ilo + DXo]} , (26) 



where 



Kr 



(27) 



The flow of k{q) follows directly from the flow equation 
for Hhh- The flow equation is shown diagrammatically in 



Fig. [2] and is given by the expression 



dA'Ehhik) = k^dAKk 



k'^{2f^o 



DXo) 



[SL(k) + S^^{k)]/2 - 5,t(fc) - S^^m ■ 



(28) 




o 





FIG. 2: (Color online) Diagrams which enter the flow of the 
self-energy E^h of the out-of-plane fluctuations. The solid dot 
above the self-energy on the l.h.s. denotes a derivative w.r.t. 
A, the other symbols are defined in the caption of Fig. [T] 



The diagrams entering the flow of Tihh{k) are 



Ghh{q)^hhhhik,~k,q,^q), 



^hhik) — 
^hh{k) — 
Shlik) = / GHh{q)Gab{q') 



Gabiq)T^hLbik,-k;q,-q) 



(29a) 
(29b) 



X r(l(fc, q; -q')4iU-fc, -q; q') , (29c) 



Shhik) = / Gatiq)GM) 



X ril(fc, q'; q)r[l{-k, g'; -q) , (29d) 



where the direction of k on the right hand side of these 
equations can be chosen arbitrarily since the diagrams 
only depend on the modulus k, and q' — k + q. The 
subscript hh indicates two external /i-legs and the super- 
script indicates the composition of the internal loop, i.e. 
Sj^f^ indicates a diagram where the loop consists of one 
/i-line. Sj^'^ differs from 5^™ in that in 5";^^ the single- 



■'hh 



^hh 



scale propagator is a u-line whereas in S^^'j^ it is a /i-line. 
With the vertices given in Eqs. (|22all22e() these diagrams 
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become 



expressions for these diagrams are 



Jq 

X [2 fig, + Ao + (mo + fiq' + K' )] , (30a) 
SLik) f q'[G±{qKD - 1) + G^^{q)]i2y^ fio + Ao) , 

Jq 

(30b) 

Sl:Uk) = k'jl I q^Ghh{q'){G^{q)fil(l-y^)(2ky + qY 



G 



II (<?) [^fi-qV{(l + ky) + \q{k + qy)] ^ | , 



(30c) 



5^,\fc) = fcVl / q^Ghh{q'){Gdqr^l{l - y'){2ky + qf 

Jq 

+ G|| ((?) [2M,y(g + ky) + \{k + qy)] ' } . 

(30d) 

Note that each of the diagrams given in Eqs. (I30all30d|) 
is to leading order quadratic in k. However, their sum 
cancels to leading order exactly the first term in Eq. ([25)) 
so that the overall leading term of the fiow of Y^hh is 
indeed quartic in k. 




O 

1/2 + 1/2 





FIG. 3: (Color online) Diagrams which determine the flow of 
the self-energy Eat of the in-plane fluctuations. The symbols 
are deflned as in Figs. [T]and[21 



GhA(Q')ritL(fc,-fc;9,-9) 

Gcd{Q)^^aLik,-k,q,-q). 



(32a) 
(32b) 



Ghh{q)Ghh{q') 



X ril(g, -g'; k)T['>,i-q, q'; -k) , (32c) 
5r(fe) = / GMGefiq') 

Jq 

X r(i(fc, q, -q')rg(-fc, -q, q') , (32d) 



where q' = k + q. After projecting on transversal and 
longitudinal components, one finds from Eq. pip and the 
form of the projected self-energies in Eqs. (j21al) and (I21b[) 



k^dAfik = k\flo + D~Xo/2)dAJl + {l/2)Sl{k) 
+ (l/2)5i(fc) - Sf\k) - S'7(fc) 

- k^dAKk , (33a) 
k^dA{2fik + \k) = k^fio + Z?Ao/2)9aJ1 + {l/2)S^{k) 

+ (l/2)5[|Xfc) - S!^'\k) - Sfik) 

— k^dA^k ■ (33b) 



The projected diagram S'^(fc), with a =± or ||, is defined 
via P;^^{k)S^^{k) = Pab{t^)Sl{k) and the other projected 
diagrams are similarly defined. With y ^ k ■ q/kq, the 
transverse projected diagrams have the form 



Finally, using the flow equation of Eafc(fe), one can de- 
termine the flow of jlk and A^. The flow of Yiab{k) has 
the form (see Fig. [3]) 

dA^abik) = [(I?Ao + 2fLk + 2fLQ)e5ab 

+ 2{\k + fLk)k''k^]dAJA + (l/2)5rb(fe) 

+ {imsUk) - s':i\k) - (fc) , (31) 

where the flrst term arises from the diagram with one 
external leg coupled to the derivative of the order param- 
eter field and S1^i^{k) stands for a diagram with external 
M-lines with flavors a, b and momentum fc, and an internal 
loop which consists of one u-line. Similarly, for S'^^{k) 
the internal loop consists of two w-lines and for Sj^l^{k) 
the internal loop consists of two /i-lines. The analytic 



srik) 



S'lik)=k^ / g2G,^,^(q)[2Mo2/' + Ao] , 

Jq 

SUk) = ^ I q^[[G^{q){D-l) + G\\{q)] 



(34a) 



D 

X (i?-l)[2/ioy'-f Ao] 
+ [G^{q){D-2 + y') + G\\(q)(l-y^)\ 

x2[M,'(l+y') + V2/']}: (34b) 



D - 1 
r,2 



q^Ghh{q)Ghh{q') 



Xf-4ik + 2yq)'{l^y'). 



(34c) 



The expression of the diagram S'""(fc) is rather long and 
can be found in the appendix. The longitudinal projec- 
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tions are 



(35a) 



Sll{k)^k' / q^{[G^{q){D-l) + G\\{q)] 

Jq 

X {D - l)[2iloy^ + ~Xo] 
+ [Gi((7)(l-y') + G||(g)z/] 

X [Xkiq + ky) + 2/ife(fc + 



(35b) 



(35c) 



and the expression for Su "{k) is given in the appendix. 



IV. RESULTS 



The set of Eqs. ([28]), ((33a|) . and ((33b| form a set 

of coupled integrodiffcrential equations which we solve 
numerically. We are mainly interested in the behavior 
of the D = 2 dimensional membrane embedded in three 
dimensional space near its critical regime. For simplic- 
ity we use as initial conditions momentum independent 
coupling constants, 



KAo = 1 and /iA„ = Aa„ = Aq , 



(36) 



The initial mean field value of the order parameter Jaq 
will be tuned to reach the critical point of the crumpling 
transition. Since our approach breaks down within the 
crumpled phase, we can approach the critical point only 
from the flat phase. We begin by discussing the proper- 
ties of the flat phase. 



A. Flat phase 

The flat phase is characterized by a finite order param- 
eter J = limA-i-o J A- Asymptotically, the out-of-plane 
fluctuations are governed by the anomalous dimension 
r]f associated with the flat phase fixed point. 



Ghhik) 



(37) 



while the in-plane fluctuations are governed by the 
anomalous dimension 77", 



Gi(9)^G||(g)^g-(2+,") 



(38) 



We can extract the flow of the anomalous dimension from 
the leading q = term of «((?), 



~(0) 



(39) 



where k 



In accordance with the deriva- 



i-A — ^^q=0- 

five expansion result jl8| and recent MC simulations of 



graphene (22| we find for D = 2 the anomalous dimen- 
sion rjf K, 0.85 which yields 77" « 0.30. These values are 
also close to the SCSA result 77/ « 0.821 [iand the MC 
results 77 « 0.750(5) [2^ and 77 « 0.81(3) [30^. A typical 
flow of the anomalous dimension far away from the crit- 
ical point is shown in the upmost curve in Fig. |4l There 
is some scatter of the numerical data for 77A since t^a has 
to be extracted from the A derivative of the very small q 
dependence of kq . The numerical noise in the function kq 
itself is however very small. The anomalous scaling of the 
propagators, as expressed through Eqs. (|57)) and is 
observable only for momenta smaller than the Ginzburg 
scale qc- This scale can be obtained perturbatively and 
has for £> = 2 the form [H, i] 



[3Xo/(27r)]i/V(4^) 



where 



Xo = 4J^A(A + A)/(2A + A) 



(40) 



(41) 



is the (bare) Young's modulus. Since for all our calcula- 
tions we used Eq. ([36]) and similar values of J^^ « 0.5, 
one flnds for all the data we present here to a good 
approximation a common Ginzburg scale of the order 
of qc w O.2A0. For q <^ qc the Green's functions 
are strongly renormalized and the flat phase fixed point 
scaling regime appears. This behavior is clearly seen 
in both Ghhil) and G±{q), see the lowest curves in 
Figs. [5] and [SI which where obtained for J^^ — 1/2 or 
SJ — Jaq — Jc ~ 0.08, where Jc is the critical Jaq value 
of the crumpling transition for the initial values of the 
elastic constants stated in Eq. ([36]) . In the small q limit, 
we find Ghh{q) oc l/q^~^f , as expected from the analysis 
of the flow of 77A . Similarly, for the in-plane correlation 
functions we find Ga oc Xjq^^"^ with a =±, || and an 
anomalous exponent 77" = 2 — 2?;/, as expected from the 
Ward identity Eq. 0. 

The upper curve in Fig. [7] shows the flow of the Pois- 
son's ratio 



A 



CTA 



g=0 



2/iq=0 + Aq=0 



(42) 



which in the flat phase is known to acquire a negative 
value for A ^ 0. For A — one finds in the flat phase 
a Poisson's ratio cr/ « —1/3, as in the NPRG deriva- 
tive expansion |18||. and in perfect agreement with both 
the SCSA result Q and MC results for phantom mem- 
branes 3f] and in good agreement with MC results for 
self-avoiding membranes [321] . 



B. Behavior close to and at the crumpling 
transition 

Upon lowering the initial value of J\^-, we can tune 
the membrane towards the critical regime of the crum- 
pling transition. Exactly at the crumpling transition the 
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0.9 




FIG. 4: Flow of the critical exponent rj for different initial 
values of J^^ (solid dots). For finite = limA-)-o 
saturates at the flat phase fixed point value 77/ ~ 0.85, this 
value is indicated by the upper vertical line. At the crumpling 
transition (with Jaq = Jc) tja saturates at the critical fixed 
point value rjc f« 0.638 (see curve 1), this value is indicated 
by the lower vertical line. The values for (Ja^ — Jc)/Jc a-re 
8.0 X iO"'^ (curve 2), 3.0 x fO"^ (curve 3), 5.0 x fO"" (curve 
4), 9.f X fO"^ (curve 5), and i.9 x fO"^ (curve 6). 

value of the anomalous dimension is changed from the 
flat phase fixed point value. The lowest curve in Fig. |3] 
shows the flow of ryA at the critical value of . One can 
clearly observe a fixed point value rjc which is different 
from the one obtained for the flat phase. We find a value 

r?c ~ 0.63(8) , (43) 

which is larger than the SCSA result rjc — 0.535 Q, 
and the value rjn = 0.47 from a MC renormalization 
group analysis , but rather close to the large D result 
77 = 2/3 and the derivative expansion approach to the 
NPRG [H] 77 w 0.627, obtained with a sharp cutoff. In 
the derivative expansion, a weak dependence of the criti- 
cal properties on the form of the regulator was reported. 
In our numerical approach, we are for reasons of numer- 
ical stability restricted to analytical regulators [s^. We 
would expect however a much weaker dependence on the 
form of the regulator in our approach since we keep the 
entire momentum dependence of all two-point vertices 
and since the regulator, for a given A, only affects the 
momentum dependence of the two-point functions. 

A small distance — > away from the critical 
value J^, rjA initially approaches the critical fixed point 
value rjc but deviates from it in the IR limit to finally 
saturate at the flat phase fixed point. This behavior can 
be seen for different values of J^^ > Jc in Fig. HI 

The behavior of the Green's function Ghh{q) at the 
critical point is shown in Fig.[5j There is again a crossover 
from perturbative scaling to anomalous scaling near qc 
similar as in the flat phase. The anomalous scaling is 



10^ 

T 
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critical FP scaling 


CO 


flat FP scaling ~ " 








pert, scaling 








(D 1 







1 0''* 10'^ 10"^,,. 1 0"^ 1 



FIG. 5: The correlation function Ghhik) (solid lines) for dif- 
ferent values of — (curves are multiplied with fc^ for 
easier visibility of the different scaling regimes) . At the criti- 
cal point of the crumbling transition (upmost curve) a single 
crossover from the perturbative scaling to a critical scaling 
with 77c ~ 0.63(8) is observed. Slightly away from critical- 
ity (middle curve), Ghh enters at small momenta the flat 
phase fixed point scaling with 77/ « 0.85. Even further away 
(lower curve), the intermediate regime, where Ghh obeys crit- 
ical scaling, disappears and a direct crossover from perturba- 
tive to fiat phase fixed point scaling is observed. Dashed lines 
indicate the different scaling behaviors. 

now the critical scaling with rj = rjc- As expected from 
an analysis of the flow of the critical exponent, there is a 
second crossover at a smaller scale qc at which the critical 
scaling regime terminates and scaling associated with the 
flat phase flxed point critical exponent 77/ appears. A 
similar two-parameter scaling is in fact also present in 
<I>^ models at weak coupling [34] . 

The scale qc can be estimated from the simple obser- 
vation that the crossover occurs for kq^q^ ~ J^jiq^. Since 
kq (X 5^'''= and flq oc q*'^^"^^", we find for D = 2 the 
relation 

qc oc J^/^^ , (44) 

where J is the fully rcnormalized magnitude of the order 
parameter which vanishes at the critical point. 

The in-plane fluctuations show a similar behavior, see 
Fig. [6] (we only show results for G±{k), but G\i (fc) shows 
a very similar behavior) . Directly at the crumpling tran- 
sition the order parameter J vanishes and one therefore 
has Ghh{q) = G±{q) = G\\{q) for all q. Similarly, a finite 
distance away from the critical surface, but for momenta 
in the regime qc ^ q ^ qc the fluctuations are still 
dominated by the vicinity of the critical fixed point and 
consequently the in-plane modes scale anomalously with 
the same exponent as Ghh, i-e., with a =_L, ||, 

Ga{q) ^ Ghh{q) ioTqc^q<^qG- (45) 
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FIG. 6: The correlation function G±{k) (solid lines) for dif- 
ferent values of J^^, — Jc (curves are multiplied with for 
easier visibility of the different scaling regimes). At the crit- 
ical point of the crumbling transition (upmost curve) a sin- 
gle crossover from the perturbative scaling to a critical scal- 
ing with rjc ~ 0.63(8) is observed. At criticality, we have 
Ghh{k) = G±{k) = G\\{k) since the order parameter J van- 
ishes. Slightly away from criticality (middle curve), Ghh en- 
ters at small momenta the flat phase fixed point scaling with 
rj^ = 4: — D — 2rif ~ 0.30. Even further away (lower curve), the 
intermediate regime, where Ghh obeys critical scaling, disap- 
pears and a direct crossover from perturbative to flat phase 
fixed point scaling is observed. Dashed lines indicate the dif- 
ferent scaling behaviors. 

On the other hand, for very small q one obtains 
G 

a 

2-1-?)" ^ 

as expected, since at small A (or small q), the NPRG 
flow is away from the crumpling transition fixed point 
and towards the flat phase flxed point. Note that the 
crossover at qc takes the in-plane modes thus directly 
from a q~'^~^'^'' scaling to a q~'^~^ scaling. Thus, a regime 
where the in-plane modes scale with an anomalous di- 
mension 4 — D — 2r]c is not present. While we do observe 
a /ig ~ Ag q^i' scaling with rf^ = A — D — 2ric in D = 2, 
this exponent is not observable in the in-plane correla- 
tion functions because ~ {A/ Aq)^''^'^^" on approach- 
ing the crumpHng transition such that the contribution 
to the self-energies via J^fiqq^ ~ J^A^g^ ~ q"^'^" (with 
A ^ q) scale with the same exponent as Kg(?*. 

Rather interesting is the behavior of the Poisson's ra- 
tio, defined in Eq. (j42p . at the crumpling transition, see 
Fig. [71 It is more than twice the value found for the flat 
phase, 

(Te« -0.71(5). (47) 

The smallest possible value of cr is —1 which is achieved 
when the ratio of the bulk modulus to the shear modu- 
lus vanishes, i.e. lim/j^o(/ife + A/j)//ife — > 0. The large 




FIG. 7: Flow of the Poisson's ratio ga- In the fiat phase, we 
find a = —1/3 whereas at the crumpling transition we find 
a much larger magnitude, a ~ —0.71(5). Both fixed point 
values for a a are indicated by dashed lines. The flow in the 
middle shows a flow of oa for a membrane close to the critical 
point, with the IR properties however dominated by the flat 
phase fixed point. 

negative value of <Tc thus implies that at criticality the 
bulk modulus is very small compared to the dominant 
shear modulus. While several other approaches (e.g. the 
SCSA 8| or the derivative expansion approach to the 
NPRG [l^) would allow to extract the value of the Pois- 
son's ratio at criticality, we are not aware of any earlier 
results. If the crumpling transition is indeed of second 
order, its Poisson's ratio would be comparable to those 
of the most strongly auxetic materials (materials with 
a negative Poisson's ratio) presently known 'ss']. The 
lack of orientational order at criticality, and the result- 
ing absence of a D = 2-dimensional plane along which 
the membrane extends, complicate however the physical 
interpretation of a large negative Poisson's ratio calcu- 
lated for D = 2. The flow of cta is shown in Fig. [7] as a 
function of the IR cutoff A. The lowest curve is the flow 
towards the fixed point of the crumpling transition. Also 
shown is the flow of the Poisson's ratio for parameters 
close to the critical ones but where the flow is ultimately 
to the value c/ = —1/3 associated with the flat phase. 

The magnitude of the order parameter J scales near 
the crumpling transition with an exponent /3, 

J ^ (<5r)^ , (48) 

where 6r = Vc — fg > measures the distance to the crit- 
ical temperature, Sr oc (Tc — T)/Tc. Mean field theory 
predicts /3 = 1/2, see Eq. (g]), while our NPRG calcula- 
tion yields a significantly different value for /?, 

/3« 0.22, (49) 

which was obtained from a fit for small values of J, see 
Fig. M Note that from Eq. (gH) and Eq. (gS]) one finds 
a simple relation between the crossover scale qc and the 
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FIG. 8: Dependence of the fully renormalized order parameter 
J as a function of &r = fc — fo. The critical exponent is found 
to be ^ f» 0.22. 

distance to the crumpling transition critical point, 

(br) . (50) 

This is just the usual hyperscaling relation /3 — v(D — 
2 + r?c)/2 for £) = 2 where v is the thermal exponent 
characterizing the divergence of the correlation length, 
with oc {6rY . Identifying ~ C"^? can deduce 

J/ « 0.69. (51) 

This value is in reasonable agreement with numerical 
works which find v « 0.85 [H and = 0.71(5) [s^ (older 
results on v can be found in Ref. [T2|), but is somewhat 
larger than the value found in the derivative expansion of 
the NPRG V « 0.52 When approaching the crum- 
pling transition from the ordered side, the in- and out-of- 
plane correlation functions do not decay algebraically for 
distances r 3> ^ with a characteristic scale ^ but merely 
decay with a different power-law. This is why we refer to 
this scale as a crossover scale rather than a correlation 
length. 

While our analysis shows that the crumpling transition 
is of second order, we based our analyis on the assump- 
tion that initially, at the UV cutoff Aq, all coupling func- 
tions are momentum independent. We ran a few tests 
to verify that our results are stable also in presence of a 
weak initial momentum dependence, but we cannot rule 
out that an unusual initial momentum dependence of the 
coupling constants could lead to an instability, e.g. at a 
finite momentum, in the renormalized model which would 
lead to a first order transition. It is also possible that in- 
cluding terms of third order in the stress tensor could 
modify the order of the transition. 

V. CONCLUSIONS 

To conclude, we have presented a thorough analysis of 
crystalline phantom membranes using a NPRG scheme 



which includes the full momentum dependence of the 
elastic coupling functions. It is a natural extension of 
the NPRG scheme based on a derivative expansion [l^ 
but yields significantly more information on the nature 
of the fluctuations. Since conflicting results on the order 
of the crumpling transition exist, it is important to in- 
clude as many correlations as possible in the ansatz for 
the effective average action. A-priori, it is difficult to de- 
cide which terms will be relevant near the transition since 
the anomalous dimension of the crumpling transition is 
very large and the crumpling transition could even be of 
weakly first order. Our nonlocal ansatz goes well beyond 
previous renormalization group treatments of crystalline 
membranes which relied on a finite number of coupling 
parameters and should thus yield more reliable results. 

In our approach we find a continuous crumpling tran- 
sition for physical membranes, i.e. D — 2 dimensional 
membranes embedded in d = 3 dimensional space, and 
we compute the associated critical exponents. We find an 
anomalous dimension rjc ~ 0.63(8) and a thermal expo- 
nent V ~ 0.69. An analysis of the scaling of the renormal- 
ized order parameter near the crumpling transition yields 
the critical exponent (5 ~ 0.22. Inside the flat phase we 
find an anomalous dimension rjf 0.85 which charac- 
terizes the asymptotic small momentum behavior of the 
out-of-plane fluctuations and an additional anomalous di- 
mension ?7" = 4 — Z) — 2?// « 0.30 which characterizes the 
asymptotic behavior of the in- plane fluctuations. We fur- 
ther analysed in detail the momentum dependence of the 
thermal fluctuations of the membrane at finite momenta. 
In both the flat phase and at the crumpling transition 
there is a crossover scale qQ which separates the anoma- 
lous scaling regime at small momenta from the perturba- 
tive regime. Near the crumpling transition, there is an 
additional crossover momentum scale qc oc \{T — Tc) /Tc\^ 
which separates an intermediate scaling regime, whose 
properties are determined by the crumpling transition 
fixed point, from the asymptotic small scaling regime 
where the flow is dominated by the flat phase fixed point. 

We further calculated the Poisson's ratio both at the 
crumpling transition and inside the flat phase. Inside 
the flat phase we recover the value af « —1/3 whereas 
we find a Poisson's ratio of much large magnitude at the 
crumpling transition, Uc ~ —0.71(5). 
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Appendix: Transverse and longitudinal projection of 
the diagram Sl^{k) 

Here we present the expression for the transverse and 
longitudinal projection of the diagram S"^{k) given in 
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Eq. (I32dp . which enters the NPRG flow equations of where we defined = fc + q. If we further introduce 
T,ab{k). While straightforward to evaluate, they have a 
complicated structure. The projections of S'^'^{k) can be 
written as: 



ST^JI / [Di^G^{q)G^{q') + D''^G\\{q)G\\{q') 

Jq 

+ i^i"Gi((7)G|| iq') + D{^G\^{q)G^{q')] 



(A.la) 

Sr = Jl I (D(^G^iq)GAq') + D^l^G\^{q)G\\{q') 



+ I?^"G^(g)G,|(g') + i?!^G||(g)G^(g') 



X = q k, Y ^q q' , Z ^ k ■ q' , 



(A.2) 



(A.lb) 

the functions D^'' with a,l3,j =||, ± can be written as 



D 



ill 



k'^Y'^ + q^Z'^ 



Y^Z' 



2XYZ 

X ^YZ 



D-2 



2Z 



2/vAfe(x + z)y ( -- + 2— 2 - — ]{q'- — )+M{D 



q q' 



k'^q^ 



2Y 



2f,,,^^XZ [1 - — ) [q^ - ^ ) [l - 2-^ 



X^ 



Y 



r,2 



Y' 



AX' 



2XYZ + X^q'^ 



D 



+ 2fikfigY{X + Z) [ -2^2 



X XY 
2 



q'q'' 



Z 



Y'Z' 



1 - 



2X,Xq,XZ [ q 

ZY^ 



2iL 

q^k? 
_X^ 

fc2 



k'^q^ 



fc2 

/c2 



4/i,,A,-XrZ ( 1 - 



y2 

2 /2 

q'^q 



ix+zr[q"-^ 



k'^q^ 



AfLq\XYZ\\ 



Z' 



k'^q 



Ailq\,X'Y\l 



q'^k^ 



ZY^X 
~ YZ^ 



(A.3a) 

X^' 



+ Ajlkjlq' 



q' 



{X + Z)\l 



2/ifcV (X + Z)Xl'(l- 



g2fc2 



fc2(7 



/2 



A2z2g2 ( 1 _ 



+ {x + zY—[ 1 



^2 \ xy2 

+ Aflkflq {X + Z) 



V2 



1 - 



q^k"^ 



ZY 



Ailq,\—^ 



1 - 



/c2 
X2 

" fc2 



A;2g 



/2 



q'^k^ 



2flg,YZ + Xg~Xq' 

. X2 
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X2 
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y2 



XYZ 

^2 
/c2 



X I- 



+ f4 

2y2 



g2q'2 



+ 2MfcA,, (X + Z) X 1 



r2 Mfc2 - 



3 



fc2 
X2 

g2 



2^,A,.X^ 1-2 



(A.3b) 

Y' 



q' 



2XY Z - 



YX 



) 



YZ 



X2 
/c2 



^2 
fc2 



(A.3c) 
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£>i-^ = Ag' ( k^y"^ + iZ^ + 1XYZ - 4- 



9' 

y \ / X'^\ (YX Y^Z 



,/2 



1- J^)+4A.M,'^i^^ (l-2^ 



1 



^2 



,y2 



+ AfilX^—\D-2 
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+ iXqflqXYZ [D-2 
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X2 
fc2 



(A.3d) 



For the terms of the longitudinal projection one finds 



Y^ik^-^)+2YZ^X-^ 
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YZ\ 
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Y^ 
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